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g ! ABSTRACT 

In order to understand the nature of magnetic reconnection in "free space" which is free 
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from any influence of external circumstances, I have studied the structure of spontaneous 
reconnection outflow using a shock tube approximation. The reconnection system of 
this case continues to expand self-similarly. This work aims 1) to solve the structure of 
reconnection outflow and 2) to clarify the determination mechanism of reconnection rate 
■ of the "self-similar evolution model" of fast reconnection. 

Many cases of reconnection in astrophysical phenomena are characterized by a huge 
dynamic range of expansion in size (~ 10 7 for typical solar flares). Although such recon- 
^ ■ nection is intrinsically time dependent, a specialized model underlying the situation has 

04 '. not been established yet. 

The theoretical contribution of this paper is in obtaining a solution for outflow structure 
which is absent in our previous papers proposing the above new model. The outflow 
has a shock tube-like structure, i.e., forward slow shock, reverse fast shock and contact 
discontinuity between them. By solving the structure in a sufficiently wide range of 
^ ■ plasma-/3: 0.001 < (3 < 100, we obtain an almost constant reconnection rate (~ 0.05: 

this value is the maximum for spontaneous reconnection and is consistent with previous 
models) and a boundary value along the edge of the outflow (good agreement with our 
+3 • simulation result) which is important to solve the inflow region. Note that everything, 

including the reconnection rate, is spontaneously determined by the reconnection system 
itself in our model. 

•i-H . 

^ ■ Subject headings: Earth — MHD — Sun: flares — ISM: magnetic fields 

1. INTRODUCTION magnetic configuration called a "current sheet 

It is widely accepted that magnetic reconnec- system". In this system, two similarly uniform 

tion very commonly plays an important role as a magnetized regions are set m contact, divided by 

powerful energy converter in astrophysical plasma a boun clary. We assume that the directions of the 

systems. However, there still remain many open magnetic field on both sides are anti-parallel to 

questions, not only regarding the microscopic c ari fy the , essence of the Problem In this case, 

physics of the anomalous resistivity, but also J he boundary carries a strong electric current, 

the macroscopic magnetohydrodynamical (MHD) hence ' we cal1 jt the current sheet - If resistivity is 

structure. In this work, our attention is fo- enhanced m this current sheet, it plays an impor- 

cused on the macroscopic evolutionary properties ^ ant role m eneT ^ conversion from the magnetic 

of magnetic reconnection, especially on the struc- form to others - In resistlve P lasm as, there are two 

ture of the reconnection outflow. fundamental processes of magnetic energy conver- 

We consider reconnection in an anti-parallel slon - 0ne 18 magnetic diffusion, and the other is 
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magnetic reconnection. 

The best known process of energy conversion in 
resistive plasmas from magnetic to other forms is 
Ohmic dissipation (magnetic diffusion). We must 
note, however, that plasmas are highly conductive 
in most astrophysical problems. Since the resis- 
tivity is very small, energy conversion by global 
magnetic diffusion takes a very long time, and 
is not applicable to many astrophysical phenom- 
ena with very violent energy releases, e.g., solar 
flares. Even in such a case, the majority believe 
that magnetic reconnection can convert the mag- 
netic energy very quickly. This is why we must 
study magnetic reconnection. 

Such energy conversion in the current sheet 
system is very important in many astrophysical 
plasma systems. The most famous example is 
the relation between solar activity and geomag- 
netospheric activity, called solar flares and geo- 
magnetospheric substorms. In this case, at least 
three current sheet systems are included. The 
first one is in the solar corona. The second and 
third are on the day-side and the night-side of 
the geomagnetosphere, respectively. Recent ob- 
servations and numerical simulations suggest that 
similar phenomena, seemingly activated by mag- 
netic reconnection, are universal, e.g., flares in ac- 
cretion disks of young stellar objects (YSOs, see 
Koyama et al. 1994, Hayashi et al. 1999), galac- 
tic ridge X-ray emissions (GRXE, see Koyama et 
al. 1986, Tanuma et al. 1999). 

We must note that many cases of actual mag- 
netic reconnection in astrophysical systems usu- 
ally grow over a huge dynamic range in their spa- 
tial dimensions. For example, the initial scale of 
the reconnection system can be defined by the ini- 
tial current sheet thickness, but this is too small 
to be observed in typical solar flares. We do not 
have any convincing estimate of the scale, but if 
we estimate it to be of the order of the ion Lar- 
mor radius, it is extremely small (~ 10° [m] in 
the solar corona). Finally, the reconnection sys- 
tem develops to a scale of the order of the initial 
curvature radius of the magnetic field lines (~ 10 7 
[m] ~ 1.5% of the solar radius for typical solar 
flares). The dynamic range of the spatial scale is 
obviously huge (~ 10 7 for solar flares). For geo- 
magnetospheric substorms, their dynamic range 
of growth is also large (~ 10 4 for substorms). 
Such a very wide dynamic range of growth sug- 
gests that the evolution of the magnetic reconnec- 
tion should be treated as a development in "free 



space", and that external circumstances do not 
affect the evolutionary process of magnetic recon- 
nection, at least at the expanding stage just after 
the onset of reconnection. 

Even if a system is not completely free from 
the influence of external circumstances, we can 
approximately treat it spontaneously evolv- 
ing system if the evolution timescale (which is 
estimated as Alfven transit time ta for the fi- 
nal system scale because the spontaneous expan- 
sion speed of the reconnection system is equiva- 
lent to the fast-mode propagation speed, see Nitta 
et al. 2001; hereafter paper 1) is much smaller 
than the timescale imposed by the external cir- 
cumstances (e.g., convection timescale). For typ- 
ical substorms, the evolution timescale of the re- 
connection is of the order of seconds or minutes 
while the convection timescale to compress the 
current sheet is of the order of hours. In such 
cases, the external circumstances simply play the 
role of triggering the onset of reconnection and 
the evolution itself is approximately free from the 
influence of the external circumstances. Though, 
of course, exceptional cases in which the external 
circumstances intrinsically influence the evolution 
can arise in particular situations (e.g., cases in 
which very fast convection drives the reconnec- 
tion), the author believes that it is worth estab- 
lishing a new model applicable to the cases which 
are free from external influences. 

However, no reconnection model evolving in a 
free space has been studied. We should note that 
most previous theoretical and numerical works on 
reconnection treated it boundary problem 

strongly influenced by external circumstances. 

In actual numerical studies, there is a serious 
and inevitable difficulty: For magnetic reconnec- 
tion to be properly studied, the thickness of the 
current sheet must be sufficiently resolved by sim- 
ulation mesh size. Usually, the thickness of the 
current sheet is much smaller than the entire sys- 
tem. Hence, in previous works, we have been 
forced to cut a finite volume out of actual large- 
scale reconnection systems for numerical studies. 

Of course, we wish to realize that the boundary 
conditions of this type of finite simulation box re- 
produce the evolution in unbounded space. How- 
ever, we should note that, in actual simulations, 
the boundary of such a simulation box necessarily 
affects the evolution inside the box, even if we use 
so-called "free" boundary conditions. This has an 
adverse influence on the evolutionary process of 
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the numerically simulated magnetic reconnection: 
When the reconnection proceeds and physical sig- 
nals propagating from the inner region cross the 
boundaries of the simulation box, the subsequent 
evolution is necessarily affected by the boundary 
conditions. This is because the information prop- 
agated through the boundary is completely lost 
and such an artificially cut-out simulation box 
will never receive the proper response of the outer 
regions. Additionally, numerical and unphysical 
signals emitted from the boundary may disturb 
the evolution. Such artificially affected evolution 
is obviously unnatural, and the resultant station- 
ary state should differ from actual reconnection 
occurring in a free space. 

Even in previous numerical studies aimed at 
clarification of the time evolution of reconnection 
(for example, series of studies originated by Ugai 
& Tsuda 1977 or Sato & Hayashi 1979), the evolu- 
tion could not be followed for a long time. This is 
mainly owing to restriction of the size of the sim- 
ulation box, hence, application of the results was 
limited to spatial scales typically, say, a hundred 
times the spatial scale of the diffusion region. We 
are interested here in an evolutionary process in 
free space without any influence of external cir- 
cumstances. In such a system, the evolution and 
resultant structure would be quite different from 
these previous numerical models. 

The same problem also appears in previous 
theoretical works. We must realize that, in 
many cases, astrophysical magnetic reconnection 
is essentially a non- stationary process because it 
grows in a huge spatial dynamic range. Most of 
these works, however, for mathematical simplic- 
ity, treat a stationary state of the reconnection in 
a finite volume (for example, see Priest & Forbes 
1986). These are obviously boundary problems 
and solutions must be influenced by boundary 
conditions. The problem is that we do not know 
how we should set the boundary condition in or- 
der to simulate external influences. In general, it 
is impossible to set. Hence, the situations argued 
in previous works are rather artificial and unnat- 
ural in many cases of astrophysical application. 

As will be discussed later (see section 5.6), if 
we cut the central region out from the entire ex- 
panding system, the central region tends to be a 
stationary state and is very similar to the well- 
known Petschek model, i.e., the system is char- 
acterized by the fast-mode rarefaction dominated 
inflow and the figure-X-shaped slow shock. The 



author thinks that stationary solutions are mean- 
ingful as the interior solution of the entire evolving 
system in free-space though such stationary mod- 
els are frequently treated as externally "driven" 
processes by boundary conditions. 

We demonstrated the above process of self- 
similar evolution model of fast reconnection by 
numerical simulation (see paper 1) and semi- 
analytic study (see Nitta et al. 2002; hereafter 
paper 2). Let us overview such an evolutionary 
process in free-space. 

We suppose a two-dimensional equilibrium state 
with an anti-parallel magnetic field distribution, 
as in the Harris solution. When magnetic diffu- 
sion takes place in the current sheet by some local- 
ized resistivity, magnetic reconnection will occur, 
and a pair of reconnection jets is ejected along the 
current sheet. This causes a decrease in total pres- 
sure near the reconnection point. Such informa- 
tion propagates outward as a rarefaction wave. In 
a low-/3 plasma (f3 <C 1 in the region very distant 
from the current sheet [asymptotic region]; as typ- 
ically encountered in astrophysical problems) , the 
propagation speed of the fast-magnetosonic wave 
is isotropic, and is much larger than that of other 
wave modes. Thus, information about the de- 
creasing total pressure propagates almost isotrop- 
ically as a fast-mode rarefaction wave (hereafter 
FRW) with a speed almost equal to the Alfven 
speed Vao in the asymptotic region. Hence, the 
wave front of the FRW (hereafter FRWF) has a 
cylindrical shape except near the point where the 
FRWF intersects with the current sheet. When 
the FRWF sufficiently expands, the initial thick- 
ness of the current sheet becomes negligible com- 
pared with the system size VAot, where t is the 
time from the onset of reconnection. In such a 
case, there is only one characteristic scale, i.e., 
the radius of the FRWF (VAot), which linearly 
increases as time proceeds. This is just the con- 
dition for self-similar growth. 

In our semi-analytic work (paper 2), the bound- 
ary condition along the edge of the outflow is 
artificially imposed approximating the result of 
our numerical simulation in paper 1. We set 
the boundary condition at y = for simplicity 
because the outflow is very thin. This bound- 
ary condition denotes the junction condition to 
the reconnection outflow. Hence, an important 
question "how is this boundary condition sponta- 
neously determined by the reconnection system" 
still remains. Our motivation of this work is to 
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clarify the mechanism to determine this bound- 
ary condition. 

The boundary condition at y = determines 
the inflow region. Hence, the reconnection rate 
which is usually denoted by the normalized in- 
flow speed toward the diffusion region is also de- 
termined. We must note that, in our self-similar 
evolution model, everything including the recon- 
nection rate is spontaneously determined by the 
reconnection system itself. We try to clarify how 
the system regulates the reconnection. 

This paper is organized as follows. We state a 
model of reconnection outflow as a kind of shock 
tube problem in section 2. Basic equations for 
outflow structure and the numerical procedure to 
solve the equations are listed in section 3. Prop- 
erties of the result, e.g., the reconnection rate 
and boundary condition along the slow shock are 
shown in section 4. In section 5, we summarize 
our study and discuss spontaneous structure for- 
mation in the reconnection outflow. 

2. MODEL OF RECONNECTION OUTFLOW 

We present a schematic picture of the reconnec- 
tion outflow in our self-similar evolution model 
(see figure 1). A similar structure of reconnection 
outflow is shown in figure 9 of Ugai 1999 as a sim- 
ulation result. Because of symmetry with respect 
to the x— and y— axes, we treat only the region 
x, y > 0. The outflow is composed of two dif- 
ferent plasmas which have different origins. The 
rear-half region (x < x c : reconnection jet) is filled 
with reconnected plasma coming from outside the 
current sheet (inflow region). The front-half re- 
gion (x > x c : plasmoid) is filled with the original 
current sheet plasma. These two regions are di- 
vided by a contact discontinuity located at x — x c . 

The entire outflow is surrounded by a slow shock 
which has a complicated 'crab-hand' shape. The 
Petschek-like figure-X shaped slow shock (oblique 
shock) is elongated from the diffusion region with 
a slight opening angle 9. There is a reverse fast 
shock (almost perpendicular shock) inside the re- 
connection jet (x — Xf). In front of the plasmoid, 
a forward V-shaped slow shock (oblique shock) 
forms. The opening angle and the locus (crossing 
point with x— axis) are and x s , respectively. 

The entire structure, including several disconti- 
nuities, is analogous to the one dimensional shock 
tube problem. The forward shock (V-shaped slow 
shock) and the reverse shock (fast shock) are 
formed by the collision of the reconnection jet and 



the original current sheet plasma, and propagate 
in both directions. Between these two shocks, the 
contact discontinuity forms. We approximate this 
reconnection outflow as a quasi 1-D problem in 
order to solve it analytically. Such an approxi- 
mation may be valid near the x— and y— axes, 
because the system is symmetric with respect to 
these axes. 

We focus our attention on the quasi one di- 
mensional problem of the figure-L shaped region 
near the x— and y— axes (We treat the region 
x, y — finite ^> D where D is the initial cur- 
rent sheet thickness. Note that this region appar- 
ently tends to coincide with the x— and y— axes 
in the self-similar stage which is a very late stage 
from the onset when we observe the evolution in 
a zoom-out coordinate [see paper 1].). Each re- 
gion between two neighboring discontinuities is 
approximated to be uniform. We note about the 
up-stream region p just above the X-shaped slow 
shock that the region between the slow shock and 
the separatrix field line (X-shaped field line reach- 
ing the X-point) is also approximated to be uni- 
form and the x— component of the velocity van- 
ishes. In this region, each reconnected field line 
has an almost straight shape and crosses the X- 
slow shock while each field line has a hyperbolic 
shape in the region above the separatrix field line. 

Region 3, between the contact discontinuity and 
the V-shaped slow shock, looks to be non uniform. 
Since we do not know how we can solve the two 
dimensional structure of this region analytically, 
we approximate this region to be uniform. The 
effect of this rough treatment is estimated in sec- 
tion 5.5. 

In the vicinity of the diffusion region of self- 
similar reconnection, the solution is quasi station- 
ary. Hence, the inflow speed toward the diffusion 
region must equal the magnetic diffusion speed. 
We here introduce the main assumption of this 
work: The magnetic diffusion speed is fully ad- 
justable with inflow speed into the diffusion re- 
gion. This is possible by tuning the thickness of 
the diffusion region if the resistivity is sufficiently 
large. Detailed discussion on the validity of this 
assumption is made in section 5.4. We note that 
the reconnection system in free space has only two 
parameters: the magnetic Reynolds number and 
the plasma-/? at the asymptotic region (see paper 
1). If the diffusion speed is adjustable, the in- 
flow speed is not restricted by the diffusion speed. 
Hence, the reconnection may not depend on the 
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magnetic Reynolds number (Yokoyama & Shibata 
1994). Thus, our attention will be focused on the 
dependence on the plasma-/?. We wish to em- 
phasize that the resistivity does not regulate the 
system in such a situation (It may play a pas- 
sive role) though electric resistivity is necessary 
in order to reconnect the magnetic field lines. We 
now wonder what mechanism does regulate the 
system. In order to answer this question, we con- 
sider the following shock tube-like model. We can 
treat the main part of such a problem by ideal 
(non-resistive) MHD. 

The quantities denoting the initial uniform 
equilibrium at the asymptotic region are gas 
pressure Po, mass density po and magnetic field 
strength B . Plasma-/? at the asymptotic region is 
defined as (3 = P / '(Bq/[2/i\) where \i is magnetic 
permeability of vacuum. In the rest of this paper, 
we use the normalization of physical quantities as 
in paper 2. We define units for each dimension 
as follows: (unit of velocity) = Vao = B / y/ppo 
(Alfven speed at the asymptotic region), (unit of 
the length) = VAot where t is the time from the 
onset of reconnection, (unit of mass density) = po, 
(unit of magnetic field) = B , (unit of pressure) = 
P/2-poVl. \ 

We set the 1-D shock tube-like problem as fol- 
lows (see figure 1). The system has 22 unknown 
quantities: P p , p p , v yp , B xp , B yp , 9, P 1 , p ± , v ± , 
Bi, x f , P 2 , p 2 , v 2 , B 2 , P 3 , P3, v y3 , B x3 , B y3 , <p 
and x s where P*, p*, v*, P* denote the pressure, 
density, velocity, and magnetic field, respectively 
(Note x c = v x3 = v 2 because no mass flux passes 
through the contact discontinuity). The suffices 
p, 1, 2 or 3 denote the region divided by the dis- 
continuities. The suffices x or y denote the vector 
components. Quantities 9, Xf, </> and x s denote 
the inclination of the X-shaped slow shock, the 
locus of the fast shock, the inclination of the V- 
shaped slow shock and the locus of the V-shaped 
slow shock (crossing point with a:— axis), respec- 
tively. 

These unknowns should be related to each other 
via conditions coming from the integrated form 
of conservation laws (i.e., the Rankine-Hugoniot 
[R-H] conditions) or other relations. The set of 
relations is listed in the next section. 

When (3 > 0.8, we must change the model be- 
cause the reverse fast shock does not form in the 
reconnection outflow. The criterion of (3 whether 
the reverse fast shock forms or not is in the range 
(3 = 0.7 — 0.8. Though we did not precisely 



estimate the critical value of (3, we can fairly 
match the solutions of these two schemes around 
(3 ~ 0.75. For the case including no fast shock, we 
remove Xf from the unknowns and set Qi = Q 2 
where Q* represents the physical quantities at re- 
gions 1 or 2. According to this alteration, the cou- 
pled equations of the unknowns are also altered, 
but we abbreviate their detailed forms because 
this alteration is trivial. 

3. BASIC EQUATIONS AND NUMERICAL 
APPROACH 

According to the above model of the reconnec- 
tion outflow, we obtain the following 22 equations 
for 22 unknown quantities (Detailed forms of each 
equation are listed in the appendix). 

1) relations between the pre-X-slow shock and the 
asymptotic region 

Frozen-in condition [eql] 
Polytropic relation [eq2] 

2) R-H jump conditions at the X-slow shock 
Pressure jump [eq3] 

Density jump [eq4] 

Velocity jump (parallel comp.) [eq5] 

Velocity jump (perpendicular comp.) [eq6] 

Magnetic field jump (parallel comp.) [eq7] 

Magnetic field jump (perpendicular comp.) [eq8] 

3) R-H jump conditions at the reverse fast shock 
Pressure jump [eq9] 

Density jump [eqlO] 
Velocity jump [eqll] 
Magnetic field jump [eql2] 

4) Magnetic flux conservation at the X-point 
[eql3] 

5) Force balance at the contact discontinuity 
[eql4] 

6) R-H jump conditions at the forward V-slow 
shock 

Pressure jump [eql 5] 

Density jump [eql 6] 

Velocity jump (parallel comp.) [eql 7] 

Velocity jump (perpendicular comp.) [eql8] 

Magnetic field jump (parallel comp.) [eql9] 

Magnetic field jump (perpendicular comp.) [eq20] 

7) Boundary Condition at the tip of the outflow 
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[eq21] 

8)Magnetic flux conservation all over the outflow 
[eq22] 

We can solve the majority of these equations by 
hand, and by substituting the solutions into other 
equations, we can reduce equations. Finally, nine 
equations 3, 4, 5, 6, 7, 11, 14, 21 and 22 remain as 
complicated nonlinear coupled equations for nine 
unknowns vi, B xp , B yp , 6, Pi, p±, Xf, <p and x s . 

We solve these coupled nine equations by an it- 
erative method (Newton- Raphson). In order to 
find the solution of nonlinear coupled equations, 
in general, a precise initial guess of the unknowns 
is required. Unfortunately, we cannot guess the 
solution so easily, so we adopt the following prim- 
itive method to find a good initial guess of the 
unknowns. 

We estimate them from the result of our nu- 
merical simulation in paper 1 for the case (3 = 0.2 
in the asymptotic region. Of course, the numeri- 
cal result is contaminated by numerical noise, and 
so we can treat our simulation result as simply a 
hint. 

Each equation can be reduced to a form (LHS 
of eq.)=0. When we substitute the initial guess 
of physical quantities from our simulation result 
for P = 0.2 into the unknown quantities, the left 
hand side does not vanish because these are not 
the true solution. By defining the residue as the 
sum of absolute values of the left-hand side of each 
equation, we can plot a graph of the residue as a 
function of nine unknown quantities. Then we im- 
prove the guess of each unknown to diminish the 
residue. This process will be iteratively contin- 
ued by hand until the residue settles to a value as 
small as possible (until the residue becomes two 
orders of magnitude smaller than residue for the 
simulation result). This is a very time-consuming 
procedure. Thus we can reach a good initial guess 
for the unknowns in order to use the Newton- 
Raphson method for (3 = 0.2. 

First, we start with the case (3 = 0.2. Once we 
find a converged solution of the Newton-Raphson 
procedure for the case f3 = 0.2, we treat it as 
the initial guess for the case, e.g., (3 = 0.11 or 
(3 = 0.21, and we have successively obtained the 
solutions for a range of 0.001 < j3 < 0.7. 

For the case (3 > 0.8, we replace the numeri- 
cal procedure with another one for the case of no 
reverse fast shock. When (3 = 0.8, we adopt the 



solution for (3 = 0.7 as the initial guess. Thus, we 
can continue the converging procedure to the case 
of very large (3 and obtain the result for f3 < 100. 

4. RESULT 

4.1. Reconnection rate 

We investigated the structure of the reconnec- 
tion outflow by using a shock tube approximation. 
The main issues of this work are 1) reconnection 
rate of the self-similar evolution model and 2) 
boundary condition along the edge of the outflow. 

Most interest is focused on the reconnection 
rate. Let us show the /^-dependence of the re- 
connection rate R. The reconnection rate is de- 
fined by the inflow magnetic flux per unit time: 
R = —VypB xp in our dimensionless form. Fig- 
ure 2 clearly shows that the reconnection rate 
is almost constant in a very wide range of (3. 
The reconnection rate varies in a dynamic range 
0.038 < R < 0.057 for 0.001 < (3 < 100. The 
maximum reconnection rate R = 0.057 is at- 
tained at f3 — 0.6. In the low (3 limit and the 
high (3 limit, R — > 0.050 and R — > 0.038, respec- 
tively. The range of reconnection rate obtained 
here is consistent with previous classical mod- 
els of fast reconnection, e.g., the Petschek model 
(10~ 2 < R < 10 _1 : See Petschek 1964, Vasyliunas 
1975). 

Though we can calculate for the range (3 < 
0.001 and (3 > 100 without any difficulty, it is 
seemingly less meaningful because the curve of R 
shown in figure 2 almost settles at each termi- 
nal value both for the low (3 limit and the high 
(3 limit. Since the asymptotic behaviors of both 
limits (3 — > and f3 — > oo are easily predicted, the 
author thinks our result is enough. 

4.2. Boundary condition along the slow shock 

In paper 2, the inflow region was discussed as a 
boundary problem by using the Grad-Shafranov 
(G-S) approach in which we solved a second order 
partial differential equation (the G-S equation) of 
elliptic type for the magnetic vector potential A[ 
of a perturbed field, where the perturbed field is 
a deviation from the original uniform magnetic 
field. In order to solve the G-S equation, we need 
a boundary condition along the slow shock (the 
locus of the boundary is approximated as y = 
when we put the boundary condition in paper 2 
because the opening angle of the X-slow shock is 
very small). We artificially imposed this bound- 
ary condition to be analogous with the result of 
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our numerical simulation in paper 1. Thus, an 
important question as to how we should impose 
the boundary condition for the slow shock bound- 
ary is still open. The answer is that this bound- 
ary condition is spontaneously determined by the 
shock tube problem discussed in this work. The 
perturbed magnetic vector potential A[ can be 
calculated along the edge of the slow shock (see 
figure 3) with respect to the above solutions. 

We can obtain the distribution of A[ along the 
x-axis (y = 0) from the solution. However, the 
boundary value we want is the distribution along 
the slow shock. Here we assume the shape of 
the slow shock as shown in figure 3. In order 
to cancel the original uniform magnetic field, we 
add the correction term —B y to the value on 
y — 0. From figure 3, the locus of the slow shock 
is denoted by y = xtan.9 for < x < x c , y = 
[(l — x s ) tan <f) — x c tan 0](x — x c )/(l — x c ) + x c tan 6 
for x c < x < 1. Thus, we obtain the distribution 
of A[ as A[ = B 1 (xf — x) + B 2 (x c — Xf) — B xtan.9 
for < x < Xf, A[ = B2(x c — x) — B^xtsnd for 
Xf < x < x c , —B y3 (x — x c ) — B {[(1 — x s ) tan</> — 
x c tan 9](x— x c )/ (1— x c )+x c tan6>} for x c < x < 1. 

We demonstrate the case (3 = 0.2 which we nu- 
merically simulated in paper 1, and compare the 
present analytic result with the previous numer- 
ical result in figure 4. The dotted line and the 
solid line show the simulation result of paper 1 
and the above result of this work, respectively. 
These two lines are qualitatively analogous, but 
quantitatively somewhat different. We will dis- 
cuss the difference in the discussion section 5.5. 

5. SUMMARY & DISCUSSION 

5.1. Summary 

We here briefly summarize the structure of the 
reconnection system obtained in our numerical 
simulation in paper 1. The feature of the inner 
region is very similar to fast reconnection regimes 
(Priest & Forbes 1986) of which the Petschek 
model is one member, i.e., we can find the fast- 
mode rarefaction dominated inflow and X-shaped 
slow shock. However, the entire structure shows 
properties of time-dependent reconnection as fol- 
lows. The reconnection outflow has a finite span 
which is increasing in proportion to the time from 
the onset. A V-shaped forward slow shock forms 
in the vicinity of the spearhead of the outflow. In 
region 3, plasmas are compressed by the recon- 
nection jet (piston effect), and squeezed plasmas 
are gushing out to make a vertical outflow. These 



are properly included in evolutionary process. 

In this paper, the structure of reconnection out- 
flow in the self-similar evolution model has been 
studied in a shock tube approximation. We have 
obtained the reconnection rate and the required 
boundary condition to solve the inflow region 
from the structure of the reconnection outflow. 

5.2. Self- Similarity of the reconnection outflow 

We have solved the structure of the reconnec- 
tion outflow for the self-similar evolution model 
of fast magnetic reconnection. The structure 
is spontaneously determined by the reconnection 
system itself as approximated by a kind of shock 
tube problem. The reconnection outflow is di- 
vided into several regions by discontinuities, e.g., 
forward shock, reverse shock, and contact discon- 
tinuity. We solved this shock tube problem by an 
iterative (Newton-Raphson) method. The mov- 
ing speeds of these discontinuities are constant in 
time in our result. Thus, as expected from the 
well-known shock tube problem, the reconnection 
outflow extends self-similarly. 

First, we argue for self-similarity of the evolu- 
tion of the reconnection system in free space. In 
the discussion of paper 2, we need an a priori 
statement to authorize a self-similar solution for 
the inflow region so that the reconnection outflow 
(in other words, the boundary condition along the 
slow shock) also evolves self-similarly. The va- 
lidity of this assumption was, however, unclear 
in previous papers. We must note that because 
we have found the self-similarly evolving solution 
for the reconnection outflow, we can ensure self- 
similar evolution of the entire reconnection sys- 
tem. 

Comparison of the outflow structure with a 
shock tube model is also argued in Abe & Hoshino 
(2001). They obtained a similar structure to our 
result, but there is a critical difference. The V- 
shaped discontinuity of the forward shock is iden- 
tified as an intermediate shock in their result in- 
stead of as a slow shock in our case. The author 
cannot explain the reason for the difference, but 
notes that their attention is focused on a much 
earlier stage than our case. The scale of the plas- 
moid is the order of the current sheet thickness in 
their case, but we are interested in a much larger 
scale structure at the self-similar stage. Actually, 
the numerical result of Ugai 1999 (see figure 9 of 
that paper) shows a similar structure to our re- 
sult. The later evolution of the simulation of Abe 



8 



Reconnection rate in the self-similar evolution model 



& Hoshino (2001) tends also to coincide with our 
result. 

5.3. (3 -dependence of the reconnection rate 

Second, our attention is focused on the /in- 
dependence of the reconnection rate R. In our 
result, R is almost constant in the range 0.001 < 
(3 < 100 (see figure 2). We could foresee this be- 
havior of the reconnection rate from the following 
intuition. 

We consider the low (3 limit in the asymp- 
totic region. When the plasma (3 at the asymp- 
totic region far outside the current sheet becomes 
smaller, the Petschek type X-slow shock must be 
stronger. This is explained as follows. The gas 
pressure in the up-stream region is very small rel- 
ative to the magnetic pressure for a \ow-(3 plasma, 
while in the down-stream region, the gas pressure 
should be very high to keep total pressure equilib- 
rium with the up-stream region. Hence, the total 
pressure equilibrium imposes a large gas pressure 
jump at the X-slow shock. This leads to a large 
jump in the mass density. 

Next we consider the high f3 limit in the asymp- 
totic region. Both sides of the Petschek type X- 
slow shock are filled with gas-pressure dominated 
plasmas in this case, and balanced with each other 
mainly by gas pressure. This results in a very 
small density jump at the slow shock (X — > 1). 
Note that this does not mean weak shock limit 
because the jump of the tangential component of 
velocity is significant {v\ ~ 1 while v xp = 0). The 
/5-dependence of the compression ratio X of the 
X-slow shock is shown in figure 5 and is consistent 
with the above discussion. 

The reconnection rate R is defined as R = 
—v yp B X p. First we discuss the low (3 limit. The 
inflow speed — v yp will increase as the slow shock 
strengthens (j3 decreases) for the following rea- 
son. The ejection speed of the reconnection jet is 
almost constant (~ Vao) while the mass density 
of the jet increases as the slow shock strength- 
ens (the compression ratio X increases as (3 de- 
creases). Thus, the inhalation of mass flux to- 
ward the slow shock is also strengthened. This 
results in an increase of — v yp (see figure 6) as f3 de- 
creases because the mass density at the pre-slow 
shock region is almost constant (though weakly 
rarefied by the fast-mode rarefaction wave). As 
the inflow speed — v yp increases, bending of the 
magnetic field line becomes remarkable, and B xp 
decreases (see figure 6). We can also understand 



the result for the high (3 limit in the same way 
(vice versa for the high (3 limit). Consequently, 
the product of — v yp and B xp , and hence the re- 
connection rate R is almost constant independent 
of f3. 

Dependence of the reconnection rate on the 
plasma- (3 is also considered in other works, e.g., 
Ugai & Kondoh 2001. Their attention is focused 
on an implicit contribution of (3 via a current- 
driven anomalous resistivity model. The aim of 
our study is intrinsically different from theirs. 
Needless to say, anomalous resistivity is very im- 
portant to confirm the validity of a localization of 
resistivity which is necessary to establish fast re- 
connection (see Yokoyama & Shibata 1994). We 
have concentrated on a macroscopic process which 
can be treated in the ideal MHD regime and paid 
scant attention to how we should introduce the 
anomalous resistivity. We should note that we 
have no definite information on how to make a rel- 
evant macroscopic (MHD) model of anomalous re- 
sistivity which should be treated in a microscopic 
regime. 

5.4. Relation with magnetic Reynolds number 

One may feel that our result is strange, because, 
while the reconnection needs electric resistivity, 
we treat only ideal MHD and do not consider the 
dependence on the magnetic Reynolds number. 
We here discuss a regulation process of the en- 
ergy conversion speed. As discussed in section 2, 
we assume the magnetic diffusion speed is fully 
adjustable with the spontaneous inhalation speed 
— v yp which is determined by the above shock tube 
problem. This assumption is valid for the case of 
sufficiently large electric resistivity (small mag- 
netic Reynolds number) in which the spontaneous 
inhalation speed is smaller than the maximum 
diffusion speed. Such a large resistivity may be 
realized by anomalous resistivity. We need an 
assumption for the anomalous resistivity because 
we do not have a convincing model for the mecha- 
nism and estimated value of anomalous resistivity 
(Coppi & Friedland 1971, Ji et al. 1998, Shino- 
hara et al. 2001). Note that the diffusion speed 
is adjustable by a change of the diffusion region 
thickness even if the resistivity is fixed to be con- 
stant, but if there is a lower limit for the thickness 
(e.g., simulation mesh size or the ion Larmor ra- 
dius), the diffusion speed is bound by an upper 
limit. 
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In paper 1, we suppose the Reynolds number to 
be 24.5. One may think that this is an incredibly 
large resistivity. We should be careful of what de- 
notes the length in the definition of the magnetic 
Reynolds number. Usually speaking, the mag- 
netic Reynolds number according to Spitzer resis- 
tivity for typical solar flares is extremely large: 

3 (L/10^m])(^/10 6 [K]) 3 / 2 (^/10- 3 [T]) 

(1) 

where L, T, B and n are the dimension of the en- 
tire system, the temperature, the magnetic field 
strength and the plasma number density, respec- 
tively. We must note that "the effective Reynolds 
number" which is based not on the dimension 
of the entire system but on that of the diffu- 
sion region (L ~ 10° [m]: the ion Larmor radius 
for solar corona) is important. Hence, the effec- 
tive Reynolds number R e for Spitzer resistivity is 
R e ~ 10 6 . If anomalous resistivity which is much 
larger than the Spitzer resistivity (typically as- 
sumed as ~ 10 5_6 x [Spitzer resistivity] in many 
numerical MHD simulations) takes place, the ef- 
fective Reynolds number is estimated as R e ~ 
10 0-1 . The magnetic Reynolds number adopted 
in paper 1 (which is the effective Reynolds num- 
ber: R e = 24.5) may be relevant for the case of 
anomalous resistivity. 

A localized resistivity is assumed as a model 
of anomalous resistivity at the center of the re- 
connection system. We put the thickness of the 
initial current sheet D as in paper 1. Here, D is 
supposed to be of the order of the ion Larmor ra- 
dius, hence, it is close to the minimum value of the 
thickness. The dimension of the resistive region is 
set larger than D (in paper 1, we set 2D) in order 
to deal with change of the diffusion region size. 
Thus, the effective thickness of the diffusion re- 
gion can vary between D and 2D. The maximum 
diffusion speed is achieved when the thickness be- 
comes D. 

The smallest effective magnetic Reynolds num- 
ber R e mr n is defined as 



V/D 



(2) 



where 77 is the magnetic diffusivity. For the value 
Remin = 24.5 in paper 1, the maximum diffu- 
sion speed into the diffusion region is determined 
by the Sweet-Parker model (Sweet 1958, Parker 



1963): v dmax ~ V A0 /y/ R emin ~ 0.2, while the 
spontaneous inhalation speed determined by the 
shock tube approximation is — v yp ~ 0.05. The 
diffusion speed does not limit the inflow speed 
because — v yp < Vd max - This is why the reconnec- 
tion rate of our model does not depend on the 
magnetic Reynolds number. 

If the Reynolds number is not so small (small 
resistivity), say, R em %n > 400, the diffusion speed 
is smaller than the spontaneous inhalation speed, 
and may limit the inflow speed. In such cases, 
the reconnection rate will be suppressed by the 
small diffusion speed and depends on the mag- 
netic Reynolds number in ways discussed in the 
previous literature. This problem is very interest- 
ing; however, detailed discussion for such cases is 
beyond the scope of this work. 

In this meaning, the value of R plotted in figure 
2 will be the maximum value for spontaneous re- 
connection as a function of (3 when the resistivity 
is sufficiently large (R emin < 400). The system 
cannot exceed this value by spontaneous inhala- 
tion of the inflow plasma without an external in- 
jection. 

A similar discussion about spontaneous recon- 
nection was made by Ugai in his series of works 
(e.g., Ugai & Tsuda 1979 and Ugai 1983). He 
obtained a result supporting Petschek's predic- 
tion, i.e., R oc (log_R e ) _1 . We here compare 
with the result of Ugai 1983. His interest is fo- 
cused on the steady state in a finite simulation 
box (several times the initial current sheet thick- 
ness). This is realized in a very late stage from 
the onset. At that stage, waves emitted from the 
central region propagate beyond the entire sim- 
ulation box (Ugai paid attention to a stage later 
than roughly 7 times the Alfven transit-time with 
respect to the entire size of the simulation re- 
gion). We should note that this is completely dif- 
ferent from our case. His resultant reconnection 
rate logarithmically varies in a range of 0.1-0.3 
for 2.5 < R e < 200, and the inflow speed varies 
in a range of 0.1-0.3 Vao while the maximum dif- 
fusion speed (defined by the minimum mesh size) 
varies between 20 > Vd max > 0.3. The author can- 
not explain why such a very large inflow speed is 
attained (e.g., inflow speed — v yp ~ 0.05 in our 
case), but we can see that their inflow speed is 
not so small compared with the maximum diffu- 
sion speed for R e > 10 2 . Thus, the inflow speed 
might be restricted by the diffusion speed and 
the reconnection rate might depend on the mag- 
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netic Reynolds number in Ugai's result. The au- 
thor thinks such a situation is unnatural for as- 
trophysical reconnection as discussed our works 
(e.g., section 1 of this paper). Since there is no 
essential difference between Ugai's numerical pro- 
cedure and our procedure, if one set a sufficiently 
wider simulation box than Ugai, he would obtain 
a similar result to ours. 

5.5. Boundary condition along the slow shock 

Finally, we discuss the boundary condition for 
the inflow region which remained as an open ques- 
tion in paper 2. We obtain the distribution of 
magnetic vector potential A[ for the perturbed 
magnetic field along the edge of the reconnection 
outflow (see figure 4). 

The result is roughly consistent with the re- 
sult of our numerical simulation, but it cannot 
explain the numerical result quantitatively. The 
major difference is in the region x > Xf. This 
is mainly due to the two-dimensional effects of 
the actual outflow. We here assume a quasi-one- 
dimensional model for the reconnection outflow 
for simplicity, but the actual outflow clearly has 
a two-dimensional structure. For example, the 
magnetic field lines in this region have a round- 
shaped upward convexity. If we make a precise 
model of the reconnection outflow and take ac- 
count of this two-dimensional shape, such a dis- 
crepancy will be reduced. Note that, in our ex- 
perience, the physical quantities in region 1 or 2 
which are important to determine the reconnec- 
tion rate are insensitive to the model of the slow 
shock shape, and so we do not need to be too con- 
cerned about this problem. Obviously, such fur- 
ther discussion cannot proceed without the help 
of a computer simulation, and is beyond the scope 
of this work. 

5.6. Relation to stationary models 

As noted in summary section 5.1, the central 
region of the entire expanding system tends to be 
a stationary state which is similar to the Petschek 



model. The author thinks that the Petschek-like 
stationary state is unique as the inner solution 
of spontaneously evolving reconnection systems 
in free-space with a locally enhanced resistivity 
(see paper 1). Even when an external circum- 
stance strongly influences the evolution, such the 
Petschek-like central structure will hold for a long 
duration of the order of a hundred times Alfven 
transit time with respect to the proper scale of the 
external circumstance (e.g., diameter of a mag- 
netic flux tube in case of solar flares). This is 
because the induced inflow speed is of the order 
of 10~ 2 V J 4 ( see paper 1), and it needs a long time 
to change the inflow region according to the in- 
fluence from the external circumstance (see paper 
!)• 

Other types of stationary states may occur in 
the following situations: 1) very fast plasma flow 
(~ Vao'- similar to or faster than the FRW prop- 
agation) is injected through the boundary, 2) 
proper scale of the external circumstance is not 
much larger than the initial current sheet thick- 
ness (e.g., in case of laboratory plasma), 3) resis- 
tive region is not so localized, etc. However, our 
self-similar evolution model may be applicable to 
many cases in astrophysics. 
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APPENDIX 

EQUATIONS FOR STRUCTURE OF THE RECONNECTION OUTFLOW 
The structure of the reconnection outflow is determined by the following equations. 

We assume that the region between the asymptotic region and the pre-shock region is filled with non- 
resistive plasma. Hence, the magnetic flux is frozen into the induced inflow. We also assume a polytropic 
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variation in the induced inflow because there is no violent process in the fast-mode rarefaction. Thus, 
we impose 



Frozen-in condition [eql] 

Bo = B^ 
Po P P 

and 

Polytropic relation [eq2] 

PoPo 7 = PpPp 1 ■ (A2) 

where 7 is the specific heat ratio. We assume 7 = 5/3 (monoatomic ideal gas). 

There are several discontinuities in the reconnection outflow, i.e., X-shaped slow shock, reverse fast 
shock, contact discontinuity, and forward V-shaped slow shock (see figure 1). We set jump conditions 
for both sides of each discontinuity: 

X-shaped slow shock Rankine-Hugoniot (R-H) jump conditions 
Pressure-jump [eq3] 

§ = 1 + -^(cos^ p ) 2 (X - 1) x - (cos 05^ + sin^ p ) 2 /2 

r p s p 

x[-2Vl x X + (cos9v yp r(X+ l)]/[((cos^ p ) 2 -XVl px f m ]} (A3) 



where 



Csp 



\/l P p/pp , 

V Apx = sin 9B xp + cos 9B yp ) 2 / (fip p ) , 
/i is the magnetic permeability of vacuum and X is the compression ratio. 

Density-jump [eq4] 



P — = X (A4) 
Pp 



Velocity (parallel) jump [eq5] 

Apx 



cos 6vi - v (cos 6v yp ) 2 - V\ 



sin Ovyp - v (cos 9v 2 yp - XV\ px ) 



(A5) 



where vo = v yp B xp /{B xp sin 6 — B yp cos9) is the shift speed of the de Hoffmann- Teller coordinate. 
Velocity (perpendicular) jump [eq6] 

' Sin ^=l (A6) 
cos 9y yp X 

Magnetic field (parallel) jump [eq7] 



singgi = [(cos^ p ) 2 - V\ px \X 

cos 9B xp + sin 6B yp (cos 9v yp ) 2 — XV\ px 



(A7) 



Magnetic field (perpendicular) jump [eq8] 

cos 9Bi 



— sin 9B xp + cos 9B yp 



(A8) 
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The compression ratio X is defined by the following equation (3rd order algebraic equation for X), 

{(— sin9B xp + cos 9 B yp ) 2 ''[(cos 9 B xp + sm9B yp ) 2 ('y — 1) p p (cos 9v yp ) 2 
+ (- sin 9B yp + cos 9 B yp ) 2 (27 P p - p p (cos9v yp ) 2 + 'jp p (cos9v yp ) 2 )]}X 3 
+{-p p (cos9v yp ) 2 [(- sm9B xp + cos9B yp ) 4 ('y - 1) 
+ (cos 9B xp + sin 9B yp ) 2 {j — 2)pp p (cos9v yp ) 2 
+ (— sin9B xp + cos 9 B yp ) 2 ((cos 9 B xp + sin 9B yp ) 2 (7 + 1) 
+2p(2jP p - p p (cos9)v yp ) 2 + -fp p (cos9v yp ) 2 )]}X 2 
+{pp p (cos9v yp [(cos9B xp + sm9B yp ) 2 ^ + 2(-$m9B xp + cos9B yp ) 2 {^ + 1) + p 
(2-fP p - p p (cos9v yp ) 2 + jp p (cos9v yp ) 2 )]}X 
+{-(7 + l)p 2 pl(cos9v yp ) 6 } = . 

This equation has three roots. We must choose a real positive root larger than unity. 
Reverse-fast shock R-H conditions 

Pressure-jump [eq9] 

f = (A9) 

where 

C / = 7 M 1 2 / (i-i/0) + (i-e?)/A + i 

with 

£/ = H + + 2 M( 2 - 7)(7 + l) 7 M 1 2 / )/(2//? 1 (2 - 7)) , 
I = 7 (1/A + 1) + (7 - l)7M 1 2 / /2 , 
A = Pi/(B 2 1 /(2p)) , 
Mi/ = (ui - XfV A )/c f , 



and 

Density-jump [eqlO] 
Velocity jump [eqll] 

Magnetic field jump [eql2] 



V A = Bo/y/pp^ . 

- = (MO) 
Pi 

vi-x f 



v 2 -x f 



£/ (All) 



f = if (A12) 
We must impose local magnetic flux conservation on both sides of region p and region 1. 
Magnetic flux conservation at X-point [eql3] 

v yp B xp + v x B x = (A13) 
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Force balance at contact discontinuity [eql4] 

* + fH + ^ (A14) 

V-shaped forward slow shock R-H conditions 
Pressure-jump [eql5] 

~ 1 (B cos 0) 2 -2V 2 0x X h + (x s sin 0) 2 (X, + 1 " 



P 3 = P 1 + -r(x a sin cj)Y(X h - 1) x 



X fe 2 [(x s sin 0) 2 - X^l/^j/xpo 



(A15) 



where c s o = \JjPq/pq, Vaox = \J {— sin 9B ) 2 / (ppo) and is the compression ratio 
Density-jump [eql6] 

Po 

Velocity (parallel) jump [eql7] 



X h (A16) 



cos 0(^3 - x a V A o) + sin 0^ 3 _ v 2 mQx - V / 



2 

AOx 



-VaqX s cos - X h V£ 0x 

where w m0 x = x s sin and Vao = B / ^fJTpo~. 
Velocity (perpendicular) jump [eql8] 

- sin <p{v x3 - x s V A q) + cos <f)v y3 _ J_ 
V A qx s sin X h 

Magnetic field (parallel) jump [eql9] 

cos 0^3 + sin0^ 3 = {v 2 m0x - Vl 0x )X h 
5 o cos0 v 2 mQx -X h Vl 0x 

Magnetic field (perpendicular) jump [eq20] 

- sin (j)B x3 + cos (j)B y3 _ ^ 
—B sin 

The compression ratio X h is a solution of the following third order algebraic equation: 



(A17) 



(A18) 



(A19) 



(A20) 



{(B sin 0) 2 [(B O sin 0) 2 (7 - l)p (x s sin 0) 2 + (B sin 0) 2 (2 7 P O - Po(x s sin 0) 2 + ~fp (x s sin 0) 2 )]}V| 
+{-po(^ sin0) 2 [( J B o sin0) 4 (7 + 1) + ( J B o cos0) 2 (7 - 2)pp (x s sin0) 2 
+ (5 o sin0) 2 ((fi o cos0) 2 (7 + 1) + (2 / u(2 7 P -po(a: s sin0) 2 + 1Po (x s sin0) 2 )))]}X 2 
+{pPo(^ sin 0) 4 ((A) cos 0) 2 7 + 2(B sin 0) 2 (7 + 1) + p(2^P - p (x s sin 0) 2 + -yp (x s sin 0) 2 ))}X A 

+{-(7 + l)/i 2 p[j(a: s sin0) 6 } = . 

We must choose a real positive root larger than unity. 

The tip of the reconnection outflow touches the FRWF (see figure 3), hence A[ = at this point. 
This condition reduces to the following equation. 
Boundary Condition at the tip of the outflow [eq21] 

—B y3 (l — x c ) + B x3 [(l — x s ) tan0 — :r c tan#] — B (l — x s ) tan0 = (A21) 
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From magnetic flux conservation, the injected magnetic flux must be redistributed in the reconnection 
jet. This leads to the following equation. 
Magnetic flux conservation [eq22] 

v yp B xp +[x f B 1 + (x c -x f )B 2 } = (A22) 
We assume the following trivial relations: 

Definition from the contact discontinuity 

x c = v 2 = v x3 

We can solve [eql] for p p , [eq2] for P p , [eq8] for B 1: [eq9] for P 2 , [eqlO] for p 2 , [eql2] for B 2) [eql3] 
for v yp , [eql5] for P 3 , [eql6] for p 3 , [eql 7] and [eql8] for v x3 and v y3 , [eql9] and [eq20] for B x3 and B y3 
by hand, then substitute them into other nine equations ([eq3], [eq4], [eq5], [eq6], [eq7], [eqll], [eql4], 
[eq21] and [eq22]) for the following unknowns: v 1 , B xp , B yp , 9, Pi, p 1 , Xf, (ft and x s . The only parameter 
included in this problem is the plasma-/? value at the asymptotic region. By using a Newton-Raphson 
routine, with an initial guess for these unknowns, we obtain converged solutions. The procedure to 
obtain the series of converged solutions is discussed in section 3 in detail. 
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Fig. 1. — Schematic figure of the reconnection outflow. Several discontinuities form in the outflow. We consider a 
quasi-one dimensional shock tube problem along the x— and y— axes. 

Fig. 2. — /3-dcpcndence of the reconnection rate. The reconnection rate is almost constant (~ 0.05) as a function of 
plasma-/3 at the asymptotic region far outside the current sheet. The obtained value of the reconnection rate is consistent 
with previous theoretical and numerical works, but we must note that this value is spontaneously determined by the 
reconnection system itself in a frame work of ideal (non-resistive) MHD. 

Fig. 3. — Model for the slow shock geometry. We adopt a simplified geometry of the "crab hand" shaped slow shock 
in order to calculate the distribution of A[ (perturbed magnetic vector potential). 

Fig. 4. — Boundary condition along the slow shock. The perturbed magnetic vector potential A[ for /3 = 0.2 is obtained 
along the slow shock based on the structure of the reconnection outflow (solid line) . This is important to determine the 
structure of the inflow region (see Nitta et al. 2002). The result of our previous simulation is also shown by dotted line. 
These two are analogous, but we can find a difference especially around the contact discontinuity (x = 0.87). This is a 
two dimensional effect of the actual outflow, and it shows the limit of our one dimensional approximation and simplified 
geometry of the slow shock. We must note that the value of the reconnection rate (crossing point with x = 0) is in good 
agreement with the simulation result. 
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Fig. 5. — Compression ratio of X-slow shock. Compression ratio X of the Petschck model-like X-slow shock near 
x = is plotted as a function of plasma-/3 at the asymptotic region. X increases as (3 decreases. 

Fig. 6. — Inflow speed and inflow magnetic field strength. In order to understand the behavior of the reconncction rate 
shown in figure 2, we plot the two factors (the inflow speed — v yp and the inflow magnetic field B xp ) of the reconnection 
rate. As (3 decreases, the X-shock strength increases (see discussion of this paper). This leads to the behavior of —v yp 
and B xp shown in this figure, and vice versa for the high (3 limit. Thus, the reconnection rate which is the product of 
these two quantities is almost constant independent of (5. 
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